# import the fdc dataset --------------------------------------------------

df_fdc <- read.csv("datasets/data_fdc.csv")

# import the map of turkey ------------------------------------------------

tr_map <- readRDS("datasets/data_map.RDS")

# plot the figure ---------------------------------------------------------

figure_3 <- tr_map +
  geom_point(data = subset(df_fdc, coffins > 0),
             aes(x = longitude, y = latitude, size = coffins),
             color = "darkgreen", shape = 1, stroke = 1.2) +
  coord_sf() +
  labs(x = NULL, y = NULL) +
  theme_minimal() +
  theme(panel.grid = element_blank(),
        axis.text = element_blank(),
        legend.position = "bottom") +
  scale_size_continuous(name = "Funerals", range = c(2, 5))

# display the figure ------------------------------------------------------

print(figure_3)

# save the figure ---------------------------------------------------------

ggsave(plot = figure_3, filename = "figures/figure_3.pdf", dpi = 1000,
       width = 9, height = 4, units = "in")